




NUS-786/2 
Volume 2 


FINAL REPORT FOR 
MODIFICATION OF CODES 
NUALGAM AND BREMRAD 


Volume 2: Code User Manual 


NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
Goddard Space Flight Center 
Glenn Dale Road 
Greenbelt, Maryland 

NASA Contract Number: NAS5-11781 


J. J. Steyn 


R. Huang 


May 1971 


NUS CORPORATION 
4 Research Place 
Rockville, Maryland 20850 


ASE FILE 
CORY 







NUS-786/2 


FINAL REPORT FOR 

MODIFICATION OF CODES NUALGAM AND BREMRAD 


Volume 2: Code User Manual 


For 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Goddard Space Flight Center 
Glenn Dale Road 
Greenbelt, Maryland 

NASA Contract Number: NAS5-11781 

By 

J. J. Steyn 
and 

R. Huang 
May 1971 


NUS CORPORATION 
4 Research Place 
Rockville, Maryland 20850 



TABLE OF CONTENTS 


Page 

1 . Introduction 1 

2. Code Operating Information 2 

2 . 1 General 2 

2.2.1 NUGAM2 Input Card Details 5 

2.2.2 NUGAM3 Input Card Details 14 

2.3 Code Output 23 

2.3.1 NUGAM2 Output 23 

2.3.2 NUGAM 3 Output 26 

Figures 1,2,3 and Table 1 30 

Appendix I Sample Input Listing for NUGAM2 35 

Appendix II Sample Input Listing for NUGAM3 37 

Appendix III Sample Output Listing for NUGAM2 39 

Appendix IV Sample Output Listing for NUGAM3 65 

Appendix V User Information for 05S 81 

Appendix VI User Information for SSALB Code 88 


ii 



1 . 


INTRODUCTION 


This present Volume II of the final report NUS-786 to the National Aeronautics 
and Space administration, Goddard Space Flight Center, was prepared by 
NUS Corporation under Contract NAS5-11781. It consists of a Code Users 
Manual for the NUGAM2 and NUGAM3 Monte Carlo codes described in 
Volume I. 

The NUGAM2 code determines gamma photon albedo and buildup distributions; 
it also determines K x-ray escape and finite source emission data. The 
NUGAM3 code determines inorganic scintillation pulse-height response 
distributions . 

General operating instructions for both codes are given in Section 2 Subsection 
2.1. Subsection 2.2 details NUGAM2 input. Subsection 2.3 details NUGAM3 
input. Subsections 2.4 and 2.5 describe the NUGAM2 and NUGAM3 sample 
output. Appendices I through IV consist of sample input and output listings. 
Appendix V consists of user information for the OSS code. Appendix VI con- 
sists of user information for the SSALB code described in Section 3.2 of 
Volume I. 
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2 . CODE OPERATING INFORMATION 


2 . 1 General 

The NUGAM2 and NUGAM3 codes are written in a generalized FORTRAN-IV 
level G, for the NASA-GSFC IBM 360/91 digital computer. It may be run 
after some modifications, on other computers with sufficient core size; i.e., 
the present versions require ~ 180K and 33 OK bytes (+ by 4 for approximate 
word size). There are no Sense Switch or special tape requirements. The 
code has been designed with a view to ease of translation for use on other 
computers. Double precision is used only in subroutines SC/RTA and SECTN 
in NUGAM3 by means of the IBM "*8" declaration statement. Removal of 
this declaration statement is sufficient for the large word size GDC- 6600 
computer. Input/output tapes are presently coded as 5 and 6, respectively, 
in both codes and in addition a punch card output on tape 7 in NUGAM3. 

Figure 1 shows the general flow for the data input cards . Input card 
details, order, formats, restrictions and location are given in Section 2.2. 
Card numbers are encircled and defined in the order in which they are read 
by the code. Input formats are standard FORTRAN-IV, as given in any FORTRAN 
manual. Sample input card listings for each code are presented in Appendices 
I and II. 

The code output is reviewed in Section 2.3. Appendices III and IV are sample 
output listings. They correspond to the sample input of Appendices land II. 

Debug type output may be obtained by input of N0PT = 1 on card 2. The user 
is cautioned with respect to profusion of output under this option, even 
though the codes exit after the number of histories = 100. 

Code running times on the IBM 360/91 are primarily a function of transport 
medium size and atomic number, and source photon energy. For example. 
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NUGAM2 executes at ~ 4,500 histories per minute for an iron transport 
medium,E =1.0 MeV, L = 3X and R= SX. NUGAM3 executes at ~ 8,000 
histories per minute for E = 0.662 or 2.754 MeV and a 3" x 3" Nal(T't) 
clad detector. 
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The NUGAM2 code consists of a main program and thirty-four subprograms . 
The subprogram namelist is as follows; 


AZI 

GTIS(2^ 

SIGPEK 

BKWARD 

GTSIGW 

SG^URCE 

CKMAT 

HEMI 

SPHERE 

CCfLLIS 

INDEX 

SPECTM 

dcqTmp 

KLEIN 

TA 

DFACT 

NUALGM 

TALLY 

DQfSE 

G^UTPUT 

TALLYX 

errGr 

PAIR 

TE 

FC2^WARD 

RNGEN 

TESTWT 

FY 

RG^T 

WHERE 

GENSIG 

SIGMAS 

XRAY 

GE<2^M 




The NUGAM3 code consists of a main program and thirty-two subprograms . 
The subprogram namelist is as follows: 


AZI 

GTISG^ 

SIGMAS 

BKWARD 

GTSIGW 

SG^RTA 

CKMAT 

INDEX 

SG^URCE 

CG^LLIS 

KLEIN 

SPECTM 

DCG^MP 

NUALGM 

TALLY 

DFACT 

G^UTPUT 

TALLYE 

errGr 

PAIR 

TESTWT 

FORWARD 

PAIRUP 

XPLG^T 

GAUSIN 

RNGEN 

VECTMX 

GENSIG 

RG^T 

ZIGZAG 

GEG^M 

SECTN 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card (Single Card) 


NRAND 

1 - 5 

15 

Initial random number (must be 
an odd number and different for 
each job submitted*) 

Card (Single Card) 



NGAMA 

1 - 5 

15 

Number of photon histories to 
be followed . 

NC^PT 

6-10 

15 

Output option signal 
= 0 No intermediate output 
+ 0 Intermediate output (code 
makes NGAMA = 100), 

NC2^R 

11 - 15 

15 

^ mput option for cards 
and ® : ^ 

= 0, ^d car^ 
©and® 

4= 0, e.g. 1, omit cards 
@,0©and© 

NQ(E 

16 - 20 

15 

An input option for card (© ; 
^ 0, omit card © 

= 0, read card , i.e. 

source photon energy 
and energy cut-off are 
the same as the pre- 
vious case . 

NXRAY 

21 - 25 

15 

An input option for card 
0, omit card WK 
= 0, read card 

IND 

26 - 30 

15 

Number of radial categorizations 


for forward and backscatter out- 
put data, i.e. RADI2 (on Card (^) 
divided into IND equal annular width 
zones 


* This requirement is dictated by the IBM-360 random number generating 
function: RNGEN. 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card (Single Card) , continued 


INE 

31 - 35 

15 

Number of equal energy intervals 
in escape spectrum, in the energy 
range 0 to ESPECT; INE ^ 10. 

IJKL 

76 - 80 

15 

Job termination signal 
^ 0,more cases to be run 
< 0,last case, eg. -1. 

Card (Single Card; see 

Figure 2 and Table 

1 ) 

NSQ(PT 

1-5 

15 

Source option signal (see text. 
Figure 2 and Table 1) 

RADIUS 

6-15 

F10.5 

Radius of the cylindrical medium (cm)*. 

W 

16 - 25 

F10.5 

Width of the rectangular medium (cm)* . 

B 

26 - 35 

F10.5 

Thickness of the rectangular 
medium* (cm). 

LENGTH 

36 - 45 

F10.5 

Length of the cylindrical or 
rectangular medium (cm) . 

RADI2 

Card 

46 - 55 
(Single Card) 

F10.5 

Largest albedo counting and 
transmitted photon counting 
radius, i.e. result categorization 
parameter (cm) . 

RADII 

1 - 10 

F10.5 

Collimation radius i.e. radius 
of photon beam at front face of 
transport medium (cm). 

DIST 

11 - 20 

F10.5 

Source-to-"transport medium" 
distance (cm) . 

THETO 

21-30 

F10.5 

Incident angle with respect 
to axis of medium (degree) . 


* "blank" according to the selected geometry (See Table 1). 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card ( 3 ^ 

(Single Card) , 

continued 


XO 

31 - 40 

F10.5 

1 Coordinates for point source 

YO 

41 - 50 

F10.5 

\ inside of the medium (cm); 

ZO 

51 - 60 

F10.5 

1 blank if option is not selected . 

Card 

(Single Card; 

input only if NC2(R = 

^ 0) 

NQ^M 

1 - 5 

15 

Number of laminae in the transport 
geometry (^ 9) . 

RA(1) 

6-10 

F5.3 

Thickness fraction of the first 
lamination* 

RA(2) 

11 - 15 

F5.3 

Thickness fraction of the second 

• 

• 

• 

« 

lamination. 

• 

RA(NO'M) 

• 

F5.3 

Thickness fraction of the 
NC*^M^^ lamination. 

Card (Single Card; input only if NG(E = 0 

) 

ESPECT 

1 - 10 

F10.5 

Source photon energy (MeV) . 

ECT 

11 - 20 

F10.5 

Lower energy cut-off for each 


photon history (MeV) . 

* thickness of lamina #1 total thickness 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card (Single Card; input only if NXRAY ^ 0) 


LX 

1 - 5 

15 

Lamination number; escape K 
x-rays are to be calculated 
for the, LX medium . 

EXl 

11 - 20 

F10.5 

Energy of the first selected K 
x-ray (MeV).* 

EX2 

21 - 30 

F10.5 

Energy of the second selected K 
x-ray (MeV) . * 

RXl 

31 - 40 

F10.5 

Intensity of the first selected K 
x-ray .** *** 

RX2 

41 - 50 

F10.5 

Intensity of the second selected 
K x-ray. ** 

Card (Single 

Card; input only 

if NC^R = 

0; see Figure 3) 

NINT 

1 - 5 

15 

Cross-section parameter; 

l^TTTVTT' 

energy intervals/group. 

ILC(W 

6-10 

15 

Croas-section parameter; 

2 ^ ^ is lowest energy bound 
of table (mQC^ units) . 

IHIGH 

11 - 15 

15 

Cross-section parameter; 
2 IHIGH jg highest energy 
bound of table (mQC^ units) . 


* if there are 3 laminae then LX = 1 is nearest to source. 

** see Table 2 

*** see Table 3 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card (^(Single Card), continued. 
NEL(l) 15 - 20 15 


Number of elements in the first 
lamination. 


NEL(2) 21 - 22 15 


Number of elements in the 
second lamination. 


NEL(N(2^M) 


Number of elements in the 
NGfM^^ lamination (NO^M ^ 9) . 


Card (^(Single Card) 
NE 1-5 

DENSTY 11 - 20 

ATOMNO 21-30 

ANDAW 31 - 40 


15 The number of energies for 

which cross-section data 
to be input. 

FI 0.5 The density of the medium 

• element for which cross- 

section data to be input (gm/cc) . 

FI 0.5 The atomic number of the 

medium element for which 
cross-section data to be input. 

FI 0.5 The atomic weight of the 

medium element for which 
cross-section data to be input. 


Note: Card types (6)and(7)are repeated for each element in the medium. 


NAME 

COLUMN 

FORMAT 

DESCRIPTION, PURPOSE 
AND USE 

Card Set 

(one or more cards) 



E(l) 

1 - 8 

E8.3 

First (lowest) energy for input of 
cross-section data (MeV). 

SIGPE(l) 

9-16 

E8.3 

Photoelectric cross-section for 
energy E(l), (barns/atom). 

SIGPP(l) 

17-24 

E8.3 

Pair- production cross-section 
for energy E(l) , (barns/atom). 

E(2) 

25 - 32 

E8.3 


SIGPE(2) 

33 - 40 

E8.3 


SIGPP(2) 

41 - 48 

E8.3 


E(3) 

49 - 56 

E8.3 


SIGPE(3) 

57 - 64 

E8.3 


SIGPP(3) 

65 - 72 

E8.3 


E(4) 

1 - 8 

E8.3 



E(NE) E8.3 Highest energy for input of 

cross-section data (MeV). 

SIGPE(NE) E8.3 Photoelectric cross-section 

for energy E(NE) (barn/atom) . 

SIGPP(NE) E8.3 Pair- production cross-section 

for energy E(NE) (barn/atom). 


10 


The remainder of this section consists of comments and additional explan- 
ation of the input just described. Only those items which it is felt require 
special treatment will be discussed. 

NRAND : 

Repeated use of the same input value will result in identical results, hence 
the instruction that arbitrary but differing values be input for a sequence of 
runs and run sets . Each value input must be an odd number. 

NSQ^PT, NGAMA: 


The cylindrical and rectangular transport medium geometry is defined in 
Figure 2 . Since it may be parallel beam with or without collimation and 
perpendicular or slant incidence; point source either outside or inside of the 
medium, fourteen options are available. The options NSC^PT= 1,2,3, . .. 

14 are shown in the Figure 2, and the various allowed input on Cards 
and (3^ are listed in Table 1 . 


It should be noted that for NSG^PT = 3, a fraction 1 - (W * RADI3**2), 

of the originated source photons may be "wasted" i.e. may not intercept the 
rectangular front face of the medium. The number of histories = NGAMA, 
should be increased accordingly. 

IND, RADI2 : 


Photons escape from the front or near faces of the geometry may be categorized 
as a function of radial distance from the axis . The number of category incre- 
ments may be specified as IND ^ 5, with 1 assumed by default. The radial 
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increments are equispaced with the largest = RADI2 ^ RADIUS . 


NQ^M, RA(I): 

The transport medium, length = LENGTH, may be made up of lamina 

each of length = (RA(I) * LENGTH), 1=1, NC2^M . For example if NG’M = 3 
and LENGTH = 10cm, RA(1), RA(2) and RA(3) = 0.2, 0.3 and 0.5, respectively, 
which corresponds to actual lengths of 2, 3 and 5 cm. 

EOT: 


Histories are terminated for photons where the energy is degraded below 
ECT. ECT should be input such that it is ^ See Figure 3 . 

NXRAY. EXl, EX2, RXl, RX2 , LX 

The code wUi calculate the escape of K x-rays, if the option NXRAY = 0. 

The intensities and energies of the selected K x-rays are input by the 
user as RXl, RX2 and EXl, EX2, respectively. The calculation is carried 
out for transport medium LX, e.g. if the medium has three laminations LX 
may be 1, 2 or 3. It is noted that the x-rays are those emitted by the user 
selected element of a compositional medium, e.g. Sm in Sm 203 . Element 
selection by the user is obtained through the order of input of cross-section 
data, e.g. in Sm 203 , the Sm x-ray is assumed if O cross-sections precede 
Sm cross-sections and the O x-ray is assumed if the cross-sections input 
order is reversed. Thus the cross-sections for the element of interest, 
in the lamination of interest, is input last. The user is cautioned that the 
energies EXl and EX2 must be within the cross-section energy range: 

^IL^W ^IHIGH^ where IL0W and IHTGH are defined in the next paragraph. 
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mm. ILQiW. IHIGH: 


The cross-section table generated by subroutine GENSIG consists of 
(IHIGH - IL0W) energy groups, each containing 2 ^^^^ sub-intervals. 

The total number of subintervals, over all groups, is equal to 1 + (IHIGH- 
IL(2(W) The energy width of each sub-interval within any given 

group is the same and equal to 2^ ^ , where N is the group 

interval number, beginning at N = 1, the lowest group. The energy bounds 
of group 1 are 2 ^^^^ and 2^^^^"^^ . The energy bounds of group N are 
2 IHIGH-I energy unit pertinent to this entire explanatory 

comment is m^c^ (= .51097 MeV) . An illustration of this comment is given 
in Figure 3 . Th6' relationship of the cross-sections generated with respect 
to ECT are also indicated. It should be noted that ECT ^ 2 ^^^^ and 
ESPECT ^ where both are tested in the same energy units . 


NEL(l), DENSTY, AT(2(MN(2(, ANDAW ; where the source medium consists of 
only a single element, e.g. Fe, the earlier descriptions are considered adequate. 
Where the source consists of such as a compound then further clarification is 
as follows; 


Cards @and must be repeated for each element in each lamination 
e.g. for Sm202, input data for Sm and for O. Continuing with Sm 203 
as the example, NEL(l) = 2, to indicate two elements (Sm and O); ATG^MNCTC, 
ANDAW are input as 62.0, 150.35 and 8.0, 16.0, respectively. Only the 
input values of DENSTY need reflect the number of atoms of Sm and O in 
Sm20s. DENSTY is determined as 


and 


\ 


DENSTYj 


Sm 


J sm .DENSTY 
ANDAWjg^^Q^ /Sm203 


2 * 150.35 
348.7 


*1.51 = 1.302 


DENSTY 


O 


3 * 16 
348.7 


*1.51 = 0.28 
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2.2.2 NUGAM3 Input Card Details 


DESCRIPTION, PURPOSE 

NAME COLUMN FORMAT AND USE 


Card (Single Card) 


NRAND 

1 - 5 

15 

Initial random number (must be 
an odd number and different for 
each job submitted*) . 

NC(CHAN 

6-10 

15 

Number of channels required in 
the Monte Carlo spectrum (^ 150) 
(photopeak mean pulse height). 

NPUNCH 

11 - 15 

15 

Option for punch card output of 
computed spectrum if # 0. 

IL(2(G 

16 - 20 

15 

XPLC^^T option = 0 linear plot 
¥ 0 log plot 

RESl 

21 - 30 

F10.5 ^ 
1 

^ Resolution constants for Gaussian 

‘ smearing function. ** 

RES 2 

31 - 40 

F10.5 - 


Card (Single 

Card) 



NGAMA 

1 - 5 

15 

Number of photon histories to be 
followed . 

NO'PT 

6-10 

15 

Output option signal: 

= 0 No intermediate output 
+ 0 Intermediate output (code 
makes NGAMA = 100). 

NC(R 

11 - 15 

15 

An input option signal for 
Cards (?), and : 

= 0, read Cards 

and ( 7 ) . 

f 0, omit Cards (s), 
and (?) . 


*This requirement is dictated by the IBM-360 random number generating 
function: RNGEN. 

**RES1 and RES2 correspond to K and n respectively, of Equation (27) of 
Volume 1 
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DESCRIPTION, PURPOSE 


NAME 

COLUMN 

FORMAT 

AND USE 

Card 2 

(Single Card) , continued 



NQ(E 

16 - 20 

15 

An input option for Card 
^ 0 omit Card @ 

= 0 read Card @ . 

IJKL 

76 - 80 

15 

Job termination signal: 

^ 0 more cases to be run 
< 0 last case . 

Card 

(Single Card) 



nsOpt 

1 - 5 

15 

Source option signal 
(See Figure 2 ) . 

RADIUS 

6-15 

F10.5 

Radius of the crystal (cm) . 

HT(2) 

16 - 25 

F10.5 

Nal-Crystal cylindrical length 
(cm) 

T 

26 - 35 

F10.5 

Side cladding thickness (cm). 

HT(1) 

36 - 45 

F10.5 

Front cladding thickness (cm) 

HT(3) 

46-55 

F10.5 

Photomultiplier thickness (cm). 

Card 

(Single Card) 



RADII 

1 - 10 

F10.5 

Collimation radius i.e. radius 
of photon beam at detector 
front face .* 

DIST 

11 - 20 

F10.5 

Source-to-NaI(T'f/)-crystal front 


face distance (cm) 


* may be used for parallel beam or point source 
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DESCRIPTION, PURPOSE 


NAME 

COLUMN 

FORMAT 

AND USE 

Card @ 

(Single Card) , continued . 


THETO 


21 - 30 

F10.5 

Incident angle with respect to 
crystal axis (degree) . 

XO 


31 - 40 

F10.5 


YO 


41 - 50 

F10.5 

f Coordinates of the point source 

1 inside of the detector.* 

ZO 


51 - 60 

F10.5 < 


Card @) (Single Card) 



ESPECT 


1 - 10 

F10.5 

Source photon energy (MeV) 

ECT 


11 - 20 

F10.5 

Lower energy cut-off each photon 
history (MeV) . 

Card(^ (Single 

Card: input only if NQ(R = 0; see Figure 3) 

NINT 


1 - 5 

15 

Cross-section parameter; 
2 NINT energy intervals /group. 

ILC(W 


6-10 

15 

Cross-section parameter; 
2 ILC/W lowest energy bound 
of table (moc^) units) . 

IHIGH 


11 - 15 

15 

Cross-section parameter; 

2 IHIGH highest energy 

bound of table (moc2 units) . 

NEL(l) 


16 - 20 

15 

Number of elements in the front 
and side cladding; = 0 if none. 

NEL(3) 


21 - 25 

15 

Number of elements in the 


photomultiplier; = 0 if none. 

*" blank" unless source is internal. 


16 



NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card (^(Single Card) 


NE 

1-5 

15 

The number of energies for 
which cross-section data to 
be input. 

DENSTY 

11 - 20 

F10.5 

The density of the medium 
element for which cross-section 
data to be input (gm/cc). 

ATQ^MNQ^ 

21 - 30 

F10.5 

The atomic number of the 
medium element for which 
cross-section data to be input. 

ANDAW 

31 - 40 

F10.5 

The atomic weight of the 
medium element for which 
cross-section data to be input. 

Note: 

Card types re^nd(^ are repeated for each element in the cladding 
and photomi^iplier; input for cladding only if HT(1) + T = 0 and input 
for photomultiplier if HT{3) = 0. 

Card Set 

(^(One or more cards) . 



E(l) 

1 - 8 

E8.3 

First (lowest) energy for input 
of cross-section data (MeV) . 

SIGPE(l) 

9-16 

E8.3 

Photoelectric cross-section for 
energy E(l), (barns/atom). 

SIGPP(l) 

17 - 24 

E8.3 

Pair- production cross-section 
for energy E(l), (barns/atom) . 

E(2) 

25 - 32 

E8.3 


SIGPE(2) 

33 - 40 

E8.3 
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NAME 


COLUMN 


FORMAT 


DESCRIPTION, PURPOSE 
AND USE 


Card Set (7) (O' 

ae or 

more cards) , 

continued . 

SIGPP(2) 

41 

- 48 

E8.3 

E(3) 

49 

- 56 

E8.3 

SIGPEO) 

57 

- 64 

E8.3 

SIGPP(3) 

65 

- 72 

E8.3 

E(4) 

1 - 

■ 8 

E8.3 


E(NE) 


E8.3 

Highest energy for input of 
cross-section data (MeV) . 

SIGPE(NE) 


E8.3 

Photoelectric cross-section 
for energy E(NE) (barn/atom) . 

SIGPP(NE) 

— 

E8.3 

Pair- production cross-section 
for energy E(NE) (barn/atom) . 
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The remainder of this section consists of comments and additional explan- 
ation of the input just described . Only those items which are not input in 
NUGAM2 or items which are felt require special treatment will be discussed. 


NGCHAN, ESPECT, RESl, RES2 : 

The code computes an energy absorption Monte Carlo histogram distribution 
with the full energy absorption line corresponding to the source photon energy 
(=ESPECT) located at the upper edge of channel N(2^CHAN . The computed 
Monte Carlo histogram is channel-by-channel Gaussian redistributed over 
the energy range zero to ESPECT, to ± 6CT. The Gaussian function standard 
deviation at each of the channel midpoints is determined as 


CTj = RESl * Ej ** RES2,(MeV) 

where 

E = channel midpoint energy (MeV) 
j = 1, 2 NQ^CHAN 

NSC^PT, RADIUS. HT(2), T, HT(1), HT(3), RADII, DI5T: 

The detector geometry is defined in Figure 2 , only options NSG'PT = 1,2,5, 

6,7/ 11/ 12 are available for code NUGAM3. The crystal cylindrical length and 

7 

radius are defined by HT{2) and RADIUS, respectively. 

The photomultiplier model is a cylinder of length and radius defined 

by HT(3) and RADIUS + T, respectively. T is the thickness of the side cladding 

encasing the crystal. The front face cladding radius is similarly defined by 
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RADIUS + T. The incident photons are collimated at the detector. The 
radius of collimation (RADII) may range from zero (axial incidence) to 
RADIUS + T. 

The source-to-crystal distance, DIST, is seen to be slightly larger than 
the source-to-detector distance by an account equal to the crystal front 
face cladding thickness HT(1). 

ECT: 

Histories are terminated for photons where the energy is degraded below 
ECT. ECT should be input such that it is See Figure 3. 

NINT. ILQ(W. IHIGH : 

The cross-section table generated by subroutine GENSIG consists of 

(IHIGH - ILGW) energy groups, each containing 2^^^^ sub-intervals. 

The total number of subintervals, over all groups, is equal to 1 + (IHIGH - 

ILGW) The energy width of each sub-interval within any given 

group is the same and equal to 2^ ^ , where N is the group 

interval number, beginning at N = 1, the lowest group. The energy bounds 

of group 1 are 2lLGWand 2lLQ(W + 1. The energy bounds of group N are 

2IHIGH-1 and 2IHIGH, The energy unit pertinent to this entire explanatory 

comment is mQC^ (= ,51097 MeV). An illustration of this comment is given 

in Figure 3. The relationship of the cross-sections generated with respect 

to ECT are also indicated. It should be noted that ECT ^ 2 ^^^^ and 
IHIGH 

ESPECT ^ 2 , where both are tested in the same energy units. 
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NE, SIGPE, SIGPP, E : 

The cross-section data required for input on card type @ may be obtained 
from the references *. The number of value-sets input from energy 

E(l) to E(NE), need only encompass the energy range to 

with spacing as per the references. The code generates its own internal 
cross-section table using logarithmic interpolation, A table of E, SIGPE, 
SIGPP is given in Appendix VI of the reference**. 


* references (13-17) as listed in report Volume I. 

** reference (4) as listed in report Volume I 
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NEL(I), DENSTY, ATQfMNGf, ANDAW : 


Where the transport medium (e.g. cladding, crystal, etc.) consists of 
more than a single element then further clarification is now given: 


Card types 6 and 7 must be repeated for each element in the 
medium e.g. for A-t.203*, input data for A-t and for O. Continuing 
with A-t-203 as the example, NEL(l) = 2, to indicate the two elements 
ATC^MNC2^'„ ANDAW are input as 13.0, 2 6.98 and 8.0, 16 . 0 respectively . 
Only the input values of DENSTY need reflect the number of atoms of 
At and O in Af/203 . DENSTY is determined as 


DENSTY 

Al 


2*anda w;a^ 

ANDAW 




*DENSTY)At202 


2*26.98 

101.96 


*4.0 =2.118 


and 


DENSTY/o 


3*16 

101.96 


*4.0 = 1.882 


NE. SIGPE, SIGPP, E: 


The cross-section data required for input on card type may be obtained 
from the references **. The number of value-sets input from energy 

E(l) to E(NE), need only encompass the energy range 2^^^^ to 2^^^^^, 
with spacing as per the references. The code generates its own internal 
cross-section table using logarithmic interpolation. A table of E, SIGPE, 
SIGPP is given in Appendix VI of the reference ***. 


* A-t203 is part of the cladding material used in the present work. 

** references (13-17) as listed in report Volume I 

*** reference (4) as listed in report Volume I 
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2 . 3 Code Output 


Throughout the discussion in this section, reference to the Sample Code 
Output listings given in Appendices III and IV are necessary and under- 
stood. Output which is adequately defined by headings or which is 
repetitive is either not discussed or mentioned only briefly. Output pages 
in the appendices are referred to by means of the encircled letters, e.g. , 
@, (b), etc . 

2.3.1 NUGAM2 Output (Appendix III) 

(a) This page consists of input data as well as escape energies and 
angles . 


@ This page consists of input photoelectric and pair production cross- 

section data as a function of energy in MeV units. Code computed 

scattering and total cross-section data are output as a function energy in 
2 

me units. The output data-set is repeated for each element in each 
lamina and the data-sets are repeated for each lamina in the transport 
medium . 

© These pages consist of code computed total cross-sections (cm 
as well as scattering and pair production cross-section fractions as a 
function of energy (MeV) for each lamina. The energy mesh is in accord 
with the input parameters illustrated in Figure 3. This output is repeated 
for each lamina. 


(g) This page illustrates the backward and forward scattering segmentation 
as already described in Section 2.3 of Volume I 


23 



@ These pages tabulate the backward (albedo) and forward scattering 
output as a function of angle in the noted solid angle. Unscattered, 
scattered and pair produced photons are categorized separately. Integral 
data is also shown. The output is repeated for each escape energy 
interval in accord with the input option. The output is for photon number, 
energy (MeV) and exposure (milliroentgen) . The output is for photons 
escaping within a radial distance, d, (counting radius) measured from 
the transport medium axis. 


@ This output is 
the user. 


(^repeated for the 'next' value of d, as chosen by 


This output consists of the energy integrated backward and forward 
escape as a function of angle in the noted solid angle. 

(H) These two pages consist of item (^in diagramatic representation. 

The two numbers in each segment are photon number (upper) and energy 
(lower). The A and B data refer to the two segments in which insufficient 
space is available for presenting the output. Items @ and @ are for the 
largest value of d, the counting radius chosen by the user. 


(j) These pages are similar to item ( e) except that they apply to the 
largest value of d, and they categorize single and double scatter escapes 
as well as photoelectric absorption, and they apply to the 'complete transport 
geometry', i.e. , front, back and sides. It is noted the selected solid angles 
are also different. 

(j) The output on this page consists largely of summarization and 
integrated data. Initial and Cut-off energies are in MeV units. The TO- 
TAL NO. OF COLLISIONS does not refer to a terminal classification and 
thus is primarily of statistical interest. The termination table consists 
of the following: 
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1 . ENERGY the number of histories terminated through scatter 

reducing the photon energy below the input Cut-off energy (ECT) 
threshold. Such terminations are considered as absorptions. 

2. WEIGHT the number of histories terminated through the 

weight being reduced to less than the termination threshold 
value (coded as 10 ^). Such terminations are regarded as 
absorptions . 

3. ESCAPE the number of histories terminated through TOTAL 

UNSCATTERED ESCAPE plus TOTAL SCATTERED ESCAPE. 

3A. TOTAL ESCAPING FRACTION escaping fraction per history, 

i.e. , item 3 divided by the total number of histories. 

4. ABSORBED (1. + 2. + 7.) the number of histories terminated 

through ENERGY plus WEIGHT plus PHOTOELECTRIC ABSORPTION, 
i.e. , the sum of items 1, 2 and 7. 

5. TOTAL UNSCATTERED ESCAPES — this is item ©repeated. 

5A. TOTAL UNSCATTERED FRACTION — item @ per history. 

6. TOTAL SCATTER ESCAPES the number of scattered escaping 

photons summed over all escape energies. 

6A. TOTAL SCATTERED ESCAPE FRACTION — item @ per history. 

7. PHOTOELECTRIC ABSORPTION the number of photoelectrically 

absorbed photons summed over all energies. 
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8. PAIR PRODUCTION PHOTONS — the total number of 0.51 

annihilation photons originating in pair production interactions. 

Items (|^ to are repeated for pair production photons under the label 
■TERMINATION PAIR PHOTONS.' Under the label 'TALLYCHECK' a value 
which should equate to the total number of photon histories, is output; 
the actual value may deviate negligibly due to 'round-off errors. ' The 
'TERMINATION ALBEDO PHOTONS' are determined as the fractions; 

(backs catter)/(incident). The 'TERMINATION BUILD-UP PHOTONS' are 
quantities 'transmitted.' The 'BUILD-UP FACTORS' are determined as 
the ratios: (scattered-plus-unscattered transmitted)/(unscattered 
transmitted) . 


2.3.2 NUGAM3 Output^ (Appendix IV) 


@ This page consists of self-explanatory input data as well as escape 
energies and angles. 


(b) This page consists of input photoelectric and pair production cross- 
section data as a function energy in MeV units. Code computed scattering 

f 2 

and total cross-section data are output as a function energy in me units. 
The output data-set is repeated for each element in the cladding, crystal 
and photomultiplier. 


These pages consist of code computed total cross-sections (cm ^) 
as well as scattering and pair production cross-section fractions as a 
function of energy (MeV) for cladding, crystal and 'photomultiplier.' 

The energy mesh is in accord with the input parameters illustrated in 
Figure 3. 
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@ These pages tabulate the detector backward and forward scattering 
output as a function of angle in the noted solid angle. Unscattered, 
scattered and pair produced photons are categorized separately. Integral 
data is also shown. The output is repeated for each escape energy 
interval in accord with the input option. 


@ This output consists of the energy integrated backward and forward 
escape as a function of angle in the noted solid angle. 


(|9 The output on this page consists largely of summarization and 
integrated data. Initial and Cut-off energies are in MeV units. The 
TOTAL NO. OF COLLISIONS does not refer to a terminal classification 
and thus is primarily of statistical interest. The termination table con- 
sists of the following: 

1. ENERGY the number of histories terminated through scatter 

reducing the photon energy below the input Cut-off energy (ECT) 
threshold. Such terminations are considered as absorptions. 

2 . WEIGHT the number of histories terminated through the 

weight being reduced to less than the termination threshold 
value (coded as 10 ^). Such terminations are regarded as 
absorptions. 

3. ESCAPE the number of histories terminated through TOTAL 

UNSCATTERED ESCAPE plus TOTAL SCATTERED ESCAPE. 

3A. TOTAL ESCAPING FRACTION escaping fraction per history, 

i.e. , item 3 divided by the total number of histories. 
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4. ABSORBED (1. + 2. + 7.) the number of histories terminated 

through ENERGY plus WEIGHT plus PHOTOELECTRIC ABSORPTION, 
i . e . , the sum of items 1 , 2 and 7 . 

5. TOTAL UNSCATTERED ESCAPES — self explanatory. 

5A. TOTAL UNSCATTERED FRACTION — item @ per history. 

6. TOTAL SCATTER ESCAPES the number of scattered escaping 

photons summed overall escape energies. 

6A. TOTAL SCATTERED ESCAPE FRACTION — item @ per history. 

7. PHOTOELECTRIC ABSORPTION the number of photoelectrically 

absorbed photons summed over all energies . 

8. PAIR PRODUCTION PHOTONS — the total number of 0.51 
annihilation photons originating in pair production interactions. 

Items to are repeated for pair production photons under the label 
■TERMINATION PAIR PHOTONS.” Under the label 'TALLY CHECK' a value 
which should equate to the total number of histories, is output; the actual 
value may deviate negligibly due to 'found-off errors. ' Backscattered 
and transmitted photons are output under the 'ALBEDO' and ' BUILD-UP' 
labels. 

@ This page lists the absorbed energy spectrum 'NUMBER' as a function 
of energy. The 'REDISTRIBUTED NO.' gives the absorbed energy spectrum 
weighted in accorded with the actual energy absorbed within adjacent 


28 


energy increment mid-points. The weighting is carried out at the end 
of each history. For example, if the absorbed energy corresponds to 
the increment upper energy edge, then 50% of the absorbed energy is 
tallied in the increment and 50% in the next higher increment. 'GAUSSIAN 
NO.' is the Gaussian distribution of the weighted distribution. The 
'NORMALIZED NO.' is the photopeak area normalized Gaussian spectrum. 


@ This output consists of photopeak (or escape peaks if the incident 
energy >1.02 MeV) energy and 'area.' The number of photons incident 
on the crystal is given; it may be less than the number of histories if * 
the detector is clad. The spectrum integral is given. The 'EFFICIENCY' 
is determined as the ratio: (number of photons 'detected')/(number of 
photons incident on the detector); photons are 'detected' after one or 
more collisions in the crystal. The photofraction is determined as the 
ratio: (photopeak area)/ (sum under the total spectrum). 


(I) This output is an 'on-line' plot of the 'REDISTRIBUTED NO.' Monte 
Carlo absorbed energy spectrum and the 'GAUSSIAN NO. ' spectrum 
described in item @ above. The photopeak is omitted from the plot in 
order to output the Compton continuum on a 'magnified' scale. The 
channel number and the plot y-axis 'paper scale' are given below the 
x-axis. The 'periods' indicate zero counts. The plot is primarily 
intended for qualitative purposes. 


Finally it is noted that the code outputs the absorbed energy and Gaussian 
spectra on punched cards for later possible user analysis. 
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FIGURES 1,2,3 
AND TABLE 1 
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FIGURE 1 


SIMPLIFIED FLOW DIAGRAM OF 
NUGAM2 AND NUGAM3 INPUT CARD SEQUENCE 









NSOPT = 1, 2 = 3, 4* 






= 11, 12 


= 13, 14* 


*NUGAM2 only 


FIGURE 2 

ALLOWED SOURCE AND TRANSPORT MEDIUM GEOMETRY 









c 

cn 
I u 
XI 0) 
3 C 
ui 0 



CO 

Oi 

!=) 

O 


pH 


33 


ENERGY CORRESPONDENCE BAR DIAGRAM 



CARDS (3A) and (3B) INPUT VALUES 



34 




APPENDIX I 

SAMPLE INPUT LISTING FOR NUGAM2 
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APPENDIX II 

SAMPLE INPUT LISTING FOR NUGAM3 
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APPENDIX III 

SAMPLE OUTPUT LISTING FOR NUGAM2 
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FACE ■ - TMRO BACK FACE 



TALUY SCATTERED ESCAPES 





NO* OP HtSTOftlCS INITIAU RNC«r,Y-- • CUT-OFP PNe«GV . TOTAL NO* OP COLLISIONS 



PHOTOELECTRIC ADSORPTIONS 0*0 



TERMINATION ALBEDO* 










- ARSO»BEO eNFHGY SPCCTRUM - 

fNDExcJi eNERcr eetji TO cet kmevi NUMaew redistributed no* gaussian no# NORHotizeb no* - 


0* 1 324E-01 
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’^*T9 A A6’-^i 
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0.1456E CO 
— 0vl-589 g -- ec ■ ■ = 
!>»172!E CC 
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0*21iaE00*— 
•3. 2*2516 ro 

o*2si6E or 
0*264eE -00- — 
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3 1 7flg - e e — 
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0.3442e--C0 — 
0*35756 CC 
J) * 3707E-00 — ^ ■ 
C.3B40E 00 

a r3 972e— 66 

0*4ICAE 00 


-0-» l 324E-9i — 
0*26406-01 
0*39726-01 - 
0*52966-01. 
0*66206-01 “ 
9*79446-01 
-0*92656-01 — 

0*10596 CO 

0*1192E CO 
0.1324E CO 
9*1456E 0*V- 
C.1589E 00 
^Vl-7 21 e 69 
0«16S4E 00 
3* 1966E- 90 - 
0*21 18E 00 
-9*22S1E rO“- 
0.23a3E CO 
-0-.-25 16 E 00 — 
0«2648E 00 
-0.278CE 09 
C *29136 CO 
■9*3045e'-<0~- 
9*3178E 00 
“0"c33 1 
0*34A2E 03 
-0-* 35 75E- -QO— 
C.3707E 00 
-0.3840E- <'0 “- 
0*39726 CO 

0*42376 CO 


-,33 . 
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APPENDIX V 


USER INFORMATION FOR OSS 


(Reproduced from ORNL-4160; 
Reference (5) of Volume I) 
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05S INPUT 


1. 

Card Input 


Card 

Format 

Description 

1 

n 

- ^ — 

20A4 

^ 15 

Alphameric identification 
, _ Desired number of colliding neutrons 


a . Incident neutron energy (eV) 

b. x-direction cosine of incident 
neutrons 

c. y-direction cosine 

d. 2 -direction cosine 

e. value of incident neutrons (1.0) 

4 1Z8 Initial random number. If the initial 

random number is negative or zero , then 
the machine stored number is used by 
the code as the initial random number. 
When running two or more cases which 
are identical except possibly for the 
desired number of colliding neutrons , 
different initial random numbers should 
be chosen to preserve statistical validity. 

5 I5,3F10,5 a. Geometry type « Negative for neutrons 

Incident on the flat surface of the 
system; zero or positive for neutrons 
incident on the lateral surface of the 
scintillator 

b. Radius (cm) of the scintillator 

c. Length (cm) of the scintillator 


82 



6 2E10.5 


d. Radius (cm) of the circle of inci- 
dence (used only if the geometry 
type is negative) 

a. Equivalent length (cm) in the scin- 
tillating material of the detector 

b. Equivalent length (cm) in the scin- 
tillating material of the light pipe 
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2 . Tape Intput 


The input consists of one binary tape (generated by code XSECT) consisting 
of a record of identification and, for each energy supergroup, a record 
made up of the mean free flight time for each energy subgroup, the non- 
absorption probability for each energy subgroup (presently ignored by 
05S) and the cumulative scattering pit»babiiity at each energy subgroup 
for each possible reaction type. 

The tape is prepared by the XSECT cross-section handling program of 05R 
as described in reference (51) in Volume I of this report. 



G^5S OUTPUT 

Printed output 

a. Title information 

b. Energy and direction cosines of incident neutrons 

c. Initial random number 

Qi vjeometry specifications — = _ 

e. Thicknesses of Cl detector and light pipe 

f. The upper and lower bounds (eV) of neutron energy range which 
may be considered 

g. The number of reaction types and the number of subgroups per 
supergroup 

h. The masses of the target nuclei 

i. The light table in order of ascending energy. (The'energies 
are in eV and the light values in cobalts) 

j . The energy of the incident neutrons (eV) 

k. The total number of Incident neutrons 

l. The number of reacting neutrons 

m. The number of reactions with hydrogen 

n. The number of inelastic carbon reactions 

12 9 

o . The number of C (n , oc) Be reactions 

12 

p. The number of C(n,n)30£ reactions 

12 12 

q. The number of C(n,p) B reactions 

12 11 

r. The number of C(n,pn) B reactions 

12 11 

s. The number of C(n,2n) C reactions 

t. The number of protons which escaped the scintillator 

u. The number of reactions per reacting neutron 

V. The zero-bias efficiency of the scintillator 
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w. The tabulation of th^ number of incident neutrons the total 
light of whose reactions fall within a given light interval. 
There is also a breakdown as to the type of reactions had by 
the neutrons in a given light range interval. 

X. The sums of the columns in the above tabulation 

y. The smoothing parameters 

z. The resultant pulse height distribution - - 

For the second and all subsequent cases in a given run items (f) through 
(i) of the above are omitted. 


86 



2 . 


a. The tabulation as described In item (w) of the above 

b. The pulse height distribution in log -log form for use in 
plotting the distribution by an x-y plotting device 
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APPENDIX VI 

USER INFORMATION FOR SSALB CODE 
(Reference Equation (32) of this Report, Volume I) 
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SSALB Input Card Details 


NAME 

COLUMN 

FORMAT 

DESCRIPTION, PURPOSE 
AND USE 

Card (T) (single card) 



NN 

1-5 

15 

The number of energies for 
which cross section data of 

the scatter medium to be 

input 

Card set @ (NN cards) i 

(Table for the code to interpolate the scatter 
photon cross-section) 

E(l) 

1-10 

F10.5 

First (lowest) energy for 
input of cross section data 
(MeV) 

CS(1) 

11-20 

F10.5 

Total cross -section for 
energy E(l) (cm ) 

Card set (5) (NUM cards) 



NUM 

1-10 

no 

Number of cases to be run 
(i.e. , number of times this 
card is input) 

N 

11-20 

no 

N-1 is the number of intervals 
in the numerical integration 
for each integral, e.g. , if 
N=51 then there are 50© 
and 50L determinations of the 
differential (a total of 2500) 

E0 

21-30 

F10.5 

Incident photon energy (MeV) 

R 

31-40 

F10.5 

Radius of the scatter medium (cm) 

EL 

41-50 

F10.5 

Thickness of the scatter 
medium (cm) 

U0 

51-60 

F10.5 

Total cross-section for 
E0 (cm~l) 

ELECT 

61-70 

F10.5 

Number of election/ cm for 


the scatter medium (unit of 
1Q24) 
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SSALB FORTRAN & SAMPLE INPUT DATA LISTING 
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